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Abstract 

We consider the Pauli-Fierz Hamiltonian with dynamical nuclei and investigate 
the transitions between the resonant electronic energy levels under the assumption 
that there are no free photons in the beginning. Coupling the limits of small fine 
structure constant and of heavy nuclei allows us to prove the validity of the Born- 
Oppenheimer approximation at leading order and to provide a simple formula for 
the rate of spontaneous decay. 
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1 Introduction 



In the quantum mechanical description of atoms and molecules one usually neglects 
the coupling to the radiation field and thus the possibility of emission or absorption 
of photons. The charged nuclei and electrons only interact via the static Coulomb 
interaction. Still the predictions for the spectra of atoms and molecules are in very good 
agreement with experimental data usually gathered through interaction with light. Also 
the predictions for the dynamical behavior of molecules agree with the motion observed 
e.g. in chemical reactions. The reason for the good agreement lies in the smallness of 
the fine structure constant a ~ that determines the strength of the coupling to the 
radiation field. 

It is by now well understood even on a mathematical level how the coupling to the 
quantized radiation field changes the spectrum of the Hamiltonian operator describing an 
atom or a static molecule, e.g. |BFSt IHHHt lAFFSt iFa ] IHaSe] . The quantum mechanical 
eigenstates become resonances, with energies close to the original eigenvalues, that decay 
nearly exponentially with a rate that can be computed perturbatively. 

The quantum mechanical understanding of the dynamics of molecules is based on the 
Born-Oppenheimer approximation. Roughly speaking one assumes that if the electrons 
are initially in a certain eigenstate relative to the nuclear positions (e.g. in the ground 
state), they will remain in the "same" eigenstate relative to the nuclear positions even 
when the latter change. The electronic state is "slaved" in this sense, but by energy 
conservation the electronic energy level serves as an effective potential for the motion 
of the nuclei. The validity of the approximation was proved in various versions |HaJoi 



MaSoil |SpTe[ IPST2I |MaSo2| . It is an adiabatic approximation relying on the fact that 



due to their large mass the nuclei move slowly compared to the lighter electrons. While 
transitions between different electronic levels (so-called non-adiabatic transitions) are 
possible even without coupling to the radiation field, the probabilities for such transitions 
are usually exponentially small in the adiabatic parameter and thus negligible. 

The content of this work is a mathematical analysis of molecular dynamics with 
the coupling to the quantized radiation field taken into account. Our first result is 
the validity of the Born-Oppenheimer approximation at leading order. This is of course 
expected, since the validity of the Born-Oppenheimer approximation has been confirmed 
experimentally in countless situations. Again the reason is the smallness of a which leads 
to small decay rates on the time scale set by the nuclear motion. This result is a rather 
straightforward consequence of combining the known quantum mechanical results on the 
Born-Oppenheimer approximation with standard time-dependent perturbation theory. 

The main mathematical and physical problem solved in this paper is the determi- 
nation of the rates of spontaneous emission for dynamical molecules, i.e. for situations 
where the nuclei undergo a nontrivial dynamics. Let us briefiy discuss an example 
where these rates are relevant. In Figure [T] some electronic energy levels for a di-atomic 
molecule are schematically plotted as a function of the nuclear separation R. The ground 
state energy actually behaves like for large R and thus leads to a rather small at- 
tractive force for separated atoms in the ground state, the so-called van-der-Waals force. 
One strategy to accelerate the production of dimers is to excite one of the atoms, so that 
the molecular system is in the first excited state that behaves like R~^ and thus leads 
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to a stronger attractive force. Once the nuclei come close, the system goes either into 
the ground state by spontaneous emission of a photon or the nuclei will only scatter and 
separate again. One is thus interested in the probability for spontaneous emission within 
a finite time interval while the nuclei are sufficiently close. However, this probability is 
not governed by a fixed decay rate since the electronic state and thus the lifetime of the 
resonance changes with the location of the nuclei. In particular no exponential decay 
law can be expected. 
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Figure 1: While moving in the electronic surface Ei, there is a configuration dependent 
probability for a transition to the ground state surface Eq through spontaneous emission 
of a photon. 

Our main result is an explicit time-dependent formula for the probability of sponta- 
neous decay of a dynamical molecule through emission of a photon, for finite times on 
the natural time scale of molecular dynamics. Since on this time scale the probability 
for spontaneous emission is quite small, it is far from straightforward to determine its 
leading order expression and to show that the remainder terms are even smaller. In 
particular we need to carefully separate the three time-scales given by the slow nuclear 
motion, the intermediate electronic motion and the fast photons. The main idea of our 
proof is the construction of subspaces that correspond to specific electronic states relative 
to the nucleonic configuration and momentum that are dressed by a cloud of virtual pho- 
tons. The restriction of the full dynamics to these subspaces is the Born-Oppenheimer 
approximation. A transition between two subspaces corresponds to a change in the 
electronic state with simultaneous emission or absorption of a free photon. 

Before we can make these ideas more precise, we have to explain the mathematical 
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model. For a molecule with / nuclei and r electrons the Hamiltonian is (in atomic units 
where h = 1 and c = 1) 



1 ' 1 

Hraol = Yl -0 5^ + «[^ee(|/) + V,^{x, y) + V^^{x)] 

acting on the state space "Hmoi := 'Hmic®'Hc\ '■= L"^ ® Ll^{M^) . Here x = (xi, . . . ,xi) 
denotes the configuration of the / nuclei, y = {yi, . . . ,yr) the configuration of the r 
electrons and a is the fine structure constant 

1 

a := — ^ . 

Att 137 

For notational simplicity we assume that all the nuclei have the same mass m^y^c and 
denote the electron mass by rriei- We also disregard spin as it would only complicate 
notation and not change the results. For the moment let Vee, and Van be the Coulomb 
potentials between electrons, electrons and nuclei, and nuclei respectively. 

Taking into account also the coupling to the quantized radiation field, the Hamilto- 
nian for the system becomes 

1 ' 11 1 11 

nuc el ^.^^ 

+ Hf + a[VM + Ven{x,y) + Vnn{x)] 

where pj^x '■— — iVa;^ and pj^y :— — iV^^.. It acts on the Hilbert space 

n := Hmol ® ^, 

where ^ is the photonic Fock space. 

The nuclear charge in multiples of the electron charge is denoted by Zj, H{ is the 
Hamiltonian of the free field and A\ is the quantized transverse vector potential in the 
Coulomb gauge, with a sharp ultraviolet cutoff A, needed to make the Hamiltonian H 
a well-defined self-adjoint operator. More explicitly 

A^(q) = -—Y / -^e,{k){e"^-''a{k,X)+e-''-''a*(k,X)), g e R^ 
(27r)2 jr^J\k\<A ^/2\k\ 

where a*{k,X) and a{k,X) are the standard creation and annihilation operators and 
{e\}x=i^2 are the photon polarization vectors. Note that we use a sharp ultraviolet 
cutoff just to simplify notation. All our proofs work without any changes for smooth 
cutoffs. Physically the cutofi^ is irrelevant for the problem at hand, as long as it is large 
compared to the energy of the emitted photons. 

As explained before, the validity of the Born-Oppenheimer approximation rests on 
the fact that the coupling to the field is small and we will be interested in the asymp- 
totics for small a. However, since also the Coulomb interaction depends on a (it is a 
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consequence of coupling to the field after all) the size of an atom or molecule as well 
as the electronic energy levels depend on a. To understand atoms and molecules by 
perturbation theory in a one thus switches to a-dependent units where the typical sizes 
and the typical energies are independent of a: introducing the Bohr-radius t] and the 
Rydberg-energy /i as 

r] = , fi = 2meia^ , 

one implements the change of units on the one-particle configuration space as 

and on the one-photon momentum space as 

: L^iRl, dk) L\Rl dk) , {U^(j)){k) = i2^/^(j){fik) . 

These transformations are canonically lifted to a unitary Ua on the full Hilbert space "H. 
Finally we fix the ultraviolet cutoff in units of Rydberg to some finite value Aq < oo, 

A = /iAo . 

A straightforward computation shows that in the new units the Hamiltonian becomes 

i=i 

+ ^(p,„-2vrVV/MA„(ay,))^ 



where we abbreviate 



e :-- 



+ Hi + VUy) + V,^{x, y) + Kn(a;) , (1) 

\ 1/2 




mm 

Note that even for the lightest nuclei e is already rather small, 

mp : proton mass . 

lOmp, 

The physical Hamiltonian ([T]) depends on the two small dimensionless parameters e and 
a. The smallness of e is the basis for the Born-Oppenheimer approximation in molecular 
dynamics and the smallness of a allows for a perturbative understanding of electronic 
resonances. 

Our aim is to construct a Born-Oppenheimer expansion for the resonances of the 
molecular system. To explain exactly what we mean by this statement we first recall 
some known results about the two limit cases which are contained in the Hamiltonian 
if*^'", the case e = 0, a 7^ (a molecule coupled to electromagnetic field with clamped 
nuclei) and the case e 7^ 0, a = (a molecule with dynamical nuclei but no coupling to 
the field). 
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1.1 Electronic resonances for fixed nuclei = 0, a > 0) 

When e = 0, the Hamiltonian depends parametrically on the nuclear configuration x 
and becomes 

/7^=°'°(x) =: Helix) + Hi + Ao), 
where the electronic Hamiltonian 

r 

H,i{x) := J2ply + Veeiy) + V,^{x,y) + Kn(x) (2) 

i=i 

is for every fixed x a self-adjoint operator on "Hei = L'^iM^y'). We assume that the 
spectrum of He\{x) is of the form 

a{H,i{x)) = {Eo{x),Ei{x), . . .} U oo), 

where Eq{x) < Ei{x) < E2{x) < . . . < T,{x) are eigenvalues of finite multiplicity below 
S(a;), possibly with an accumulation point at S(a;) and absolutely continuous spectrum 
in [S(a:), +00). As shown in [Zi], if 

I 

i=i 

then Hei{x) has an infinite number of eigenvalues below the threshold 

Under the same hypothesis ^ it was shown in |LiLoj . using a binding condition 
introduced in |GLLj . that H^=^''^ has a ground state E{x) for every a > (using a 
smooth ultraviolet cutoff). The existence of the ground state for small values of the fine 
structure constant a has been shown before in |BFSj . 

It is expected that the electronic eigenvalues £'i(x), £"2(2;), • • • turn into resonances 
and that apart from the ground state the spectrum of H^^^'"' is absolutely continuous. 
It was shown in |BFSj (see also [AFFSj IHHHj ) for the case I = 1 (an atom) that 
the eigenvalues {Ej}j^o become resonances in the sense of the Aguilar-Balslev- Combes- 
Simon theory ( |HiSi] chapters 16-18, |ReSi4| sections XII. 6, XIII. 10, [Sij). We quote a 



typical result (cf. e.g. Corollary 2 from |HHH] ) on the exponential decay of resonant 
states without giving technical details. 

Almost exponential decay of atomic resonances. Let Pj be the spectral projection 
of Hei = Hci{0) corresponding to the eigenvalue Ej, j 7^ 0, and let Qo be the projection 
on the Fock-vacuum. For \E' G RanPj ® Qq normalized to one it holds that 

II (P, ® Qo) e-'*^^^°'>||2 = e-*"'^ + 6(«, t) , 

where b(a,t) = 0{a^) uniformly in time and 7 > 0. It follows that the lifetime of the 
resonance is of order a^^. 

The difficult part in proving such results is to control the error term uniformly in time. 
For short times t ^ a^^ the decay rate into any other state Ei < Ej can be easily 
computed by a perturbative argument known as Fermi's golden rule. 
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Fermi's golden rule. Let Pj and Pi be the spectral projections of H^i = -ffei(O) corre- 
sponding to the eigenvalues Ei < Ej and let Qq he the projection on the Fock-vacuum. 
Then for \1/ G RanPj (g) Qo normalized to one it holds that 

uniformly on hounded time intervals, where \Dij\ is the dipole-matrix element, c.f. Sec- 
tionlTR 

Since the natural time scale for nuclear dynamics is short in this sense, we will not be 
interested in results on exponential decay on long time scales for dynamical molecules 
(it is not even clear what this would exactly mean), but in explicit decay rates in the 



form of Fermi's golden rule. However, as will be explained in Section L3, for moving 
nuclei the decay rate depends on the configuration of the nuclei, which in turn changes 
quickly on the time scale of the resonance. As a consequence one can not just adapt the 
usual perturbative argument in order to compute decay rates for dynamical molecules. 
But before coming to the full problem, let us first recall some basic facts about the case 
a = and e > 0. 



1.2 Dynamical nuclei without coupling to the field 

(£ > 0, a = 0) 

This case is the setting of the standard time-dependent Born-Oppenheimer approxima- 
tion (e.g. [HaJoil jMaSoil SpTe |Teu2[ |PSTi[ |MaSo2| ). The Hamiltonian has the form 

where we omit the field Hamiltonian Hf because it commutes with the rest and is 
therefore irrelevant. 

For kinetic energies of order one (in units of Rydberg!) the nuclei have velocities of 
order e. The time scale on which the nuclei move distances of order one (in units of 
Bohr radii) are thus times of order e~^. Hence it is natural to change the unit of time 
as well and to solve the time-dependent Schrodinger equation 



dt 

The long time-scale will be reflected in the following by the fact, that we evaluate unitary 
groups at times -, i.e. we consider e"'^^ for t of order one. 

To avoid additional technicalities, one assumes that H^i is in a suitable sense a smooth 
function of x. This requires to introduce a smearing of the nuclear charge distribution 
if G C(f^(]R^) with if > and J (f = 1. The electronic repulsion remains unchanged 



I Vn Urn \ 



n=l m=n+l 

while the electron-nucleon attraction and the nuclear repulsion become 

n=l m=l 



r I 
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l-l I „ 



and 

n=l m=n+l 

In |MaSo2| a "twisted" pseudo-differential calculus is introduced, which, generaliz- 
ing Hunziker's distortion analyticity method jHuj . allows to treat also the case of the 
unsmeared Coulomb potential. 

Let Ej{x) be an eigenvalue of the electronic Hamiltonian Hei{x) which is globally 
isolated by a gap from the rest of the spectrum. 

Definition of isolated electronic eigenvalues. Let for all x E M^' be Ej{x) an 
eigenvalue of the electronic Hamiltonian Hc\{x). The family Ej{x) is called isolated, if 
there exist two functions f± G Cb(Mi^', M) defining an interval I{x) = [/-(x), f+{x)] such 
that 

a{H^i{x)) n I{x) = Ej{x) and inf dist(Ej(x), a{H^i{x)) \ Ej{x)) > . 



This condition implies that Ej{x) and the spectral projection Pj{x) onto the eigen- 
space of Ej(x) are smooth functions of x, c.f. Lemma [4j We denote by Pj the direct 
integral 



dx Pj{x) 



which acts on n^nc ® "Hd = L\W\ Uei). 

The Born-Oppenheimer approximation rests on the observation that the electronic 
state adjusts adiabatically to the slow motion of the nuclei, i.e., that the subspace PjT-Lmoi 
is approxiamtely invariant under the time evolution. 

A rigorous version of this statement is the following theorem from |SpTe| , which is 
also a special case of Proposition [T] proven below. 

Leading order Born-Oppenheimer approximation. The operator 



:i - Pi)H 



e,a=0( 



P 



is self-adjoint on the domain D of H^'" ° and satisfies 



^(0,Wmol) 



0{e\t\ +e) . 



(5) 



This result is optimal in the sense that the difference is not smaller than order e. How- 
ever, the overlap 

||P,e-'i^^'"^"P,||=0(e) for^^j 

that the true time evolution introduces between the different electronic subspaces does 
not correspond to actual transitions between electronic states. Indeed, the subspaces 
PjUmoi can be replaced by slightly deformed superadiabatic subspaces PJ'Hmoi that are 
invariant to higher order in e. Physically in PJ'Hmoi the electronic state now depends 
also on the velocity of the nuclei. 
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We construct such superadiabatic projections in Proposition |2j A straightforward 
consequence is the following statement. 



Second order Born-Oppenheimer approximation. The operator 

Te,a=0 Tye , f -i D£\TJS,a=Oi 



+ (1 - p'^w 



is self-adjoint on the domain D of H^''^ ° and satisfies 

e ^ ^ ] l(_oo,s)(-n 



£{Wmol) 



0{e'\t\). 



(6) 



Here an energy cutoff at an arbitrary but fixed energy E is needed. This improved 
approximation of the dynamics is necessary for obtaining error terms smaller than the 
effect we are interested in, namely transitions between different electronic levels due to 
spontaneous emission of photons. But it turns out that a rigorous control of these error 
terms requires to prove (|5| and ^ with respect to more general energy norms, which 
is the main new content of Propositions [T] and |2j 

On the other hand, ^ can be shown with an error of order for any G N. 
Martinez and Sordoni even prove exponential error bounds without assuming a regular- 
ization on the nuclear charges, |MaSo2| . However, the task of computing the exponen- 
tially small transition probabilities between superadiabatic subspaces (transitions that 
happen without emission of photons) is extremely difficult even on a heuristic level, see 
|HaJoo[ [BCmiBiG^ . 

Another question is, whether one can dispose with the gap condition. At cross- 
ings of electronic eigenvalues the Born-Oppenheimer approximation breaks down and 
transitions between the levels occur at order cf. |LaTej and references therein. For 
eigenvalues embedded into or at the threshold to continuous spectrum the rate of tran- 
sition depends on the details of the model (see e.g. [Ten] and |TeTe] ) . 

Finally we remark that the importance of the Born-Oppenheimer approximation lies 
in the observation that the diagonal Hamiltonian when acting on states in the range 
of Pj, has an asymptotic expansion starting with very simple terms, 

H^^P, = {-{eVT"'? + E,{x)) P, + 0{e') . 

Here V;^*^"^ := PjV^Pj is the so-called Berry connection. Note that the electronic 
eigenvalue Ej{x) appears as an effective potential for the motion of the nuclei. To get 
the correct higher order terms, one needs to expand Hj on the range of PJ instead. 
While at zeroth and first order one obtains the same expansion as for HjPj, starting 
at second order additional terms appear, see for example in IPST2I . Thus the unitary 



groups e~^^^i Pj resp. e~^^^^ PJ can be computed by solving a Schrodinger equation for 
the nuclei only. Note that ([s]) and ^ are indeed the nontrivial mathematical statements 
to prove for justifying the time-dependent Born-Oppenheimer approximation. 



9 



1.3 Dynamical nuclei with coupling to the field 

> 0) 

The coupling to the quantized radiation field presumably turns all electronic eigenvalues 
except for the ground state into resonances. Our first aim is to prove that the Born- 
Oppenheimer approximation for a molecule described by if^'" remains valid. This makes 
sense only if the lifetime of the resonance, given according to the above discussion by 
a;~^, is bigger than the time scale of molecular dynamics, given by e^^. To control 
the relation between the two scales we thus choose a to be a function of e such that 
a{e)~^ > e~^. Assuming that 

a{e) =e^, /3 > 0, (7) 
this condition implies that /3 > |- This is always true for realistic nuclei because 

1 1 

rrip < mnuc < 250mp corresponds to Emm = ^ < £ < = e^max , 

DoU 4o 

where rrip is the proton mass. Thus 

In a 1 o In a 

/Smin = j ~ 0.75 and /3max = :j ~ 1-31 , 

which suggests to consider | < /3 < |. For some results we are able to cover even the 
range | < /3 < |, while for others we have to restrict to | < /3 < | which corresponds 
to trip < mnuc < 72mp. 

Inserting ([T]) into H^'"" as given in ([T| and expanding in powers of e, we get a 
Hamiltonian which depends just on e. Setting 

I 

m := -An'/'j2^(''y^)-p^^y 

i=i 

I 

i=i 

3=1 3=1 

where we normal ordered the quadratic terms, we can write := as 

= HI + e^^Hl + e^^+^Hl . (8) 

Note that we think of spj^x being of order e^, since we want to look at states with nuclear 
kinetic energy of order . The leading order term H^ contains no coupling to the field 
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at all. The first order term Hf describes the linear coupling of the electrons to the 
field and will be the relevant term for understanding spontaneous emission of photons. 
Contributions from i7| will always be of lower order and contribute only to our error 
terms. 

Lemma[T]below asserts that is a well-defined self-adjoint operator for e sufficiently 
small and that the expansion ([s]) makes actually sense, since the coefficients Hf and iff 
are relatively ifo'bounded with relative bounds independent of e. 

We come now to an informal statement of our main results. Let 

^^fieid := ^1 ® 1 + l®^f- 

In Corollary [T] we show that, up to a worse error estimate, the statement of ^ remains 
valid. 

Leading order BO-approximation with coupling to the field. 

For I < (3 < ^ it holds that 

= 0{e'^^-'\t\ + e) . 

c(Do,n) 

Technically this is a straightforward perturbative consequence of (|5|, as the contribution 
of e^/^Hf is of order £2^"^ for times of order e^^. However, we still believe that this result 
is conceptually important. It shows that, in the context where the Born-Oppenheimer 
approximation is usually applied, the coupling to the radiation field is negligible at 
leading order. To our knowledge this is the first mathematical result of this type. 

Our main result concerns the failure of the Born-Oppenheimer approximation be- 
cause of spontaneous emission of photons. However, as we will show, the probability for 
making a transition through spontaneous emission is of order e^'^'^ltl, which for /3 > 1 
is smaller than e"^, the square of the error in the standard Born-Oppenheimer approxi- 
mation ^ . Hence for /3 > 1 we need to consider transitions between the superadiabatic 
subspaces PJ in order to correctly separate transitions through spontaneous emission 
from errors in the adiabatic approximation. Our main result is then the following, cf. 
Theorem [2l 

Probability for spontaneous emission. Let Ej{x) > Ei{x) for all x and let \E' = 

ip with ip e RanPJ and f2 G the vacuum state. The probability for ending up in 
the i-th electronic state after time t when starting in \E' is 

\m ® l)e-'i^^ n' = e'^-' f i lllAilAf e-'t^^P,-V^||'ds + o{e'P~^) . (9) 

Here Ae and Dij are real-valued multiplication operators, namely A e{x) = Ej{x)—Ei{x) 
the difference in energy and Dij{x) := Y7k=i{Vi{^)\yk'^j{^))'Hci the dipole coupling ele- 
ment. ipi{x) and (Pj{x) are normalized electronic states in RanPj(a;) and RanPj(a;) 
respectively. 

Thus the decay probability can be computed by propagating the initial molecular 



£ — e ^ j^ficid 
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wave function according to the standard Born-Oppenheimer approximation in the level 
Ej and integrating the decay rate along this trajectory. The decay rate is given by 
I a^e~^\Dij{x)\'^AE{xY as a function of the nuclear configuration x. Recall that a = e^^ 
and that we rescaled time by e"^. This is the natural generalization of Fermi's golden 
rule for atoms, c.f. Q, to moving nuclei. 

We briefly sketch the strategy of our proof and comment on some difficulties. The 
basic idea is to use time-dependent perturbation theory according to the splitting 

= HI,,,, + el^Hl + 00^-^') . 

Note that the 0{e'^) term contains the iff term, i.e. higher order terms in the coupling to 
the field, and error terms from the Born-Oppenheimer approximation, i.e. P^H^{1 — PJ) 
and its adjoint. Abbreviating P^^ := ^ Qq, first order time-dependent perturbation 
theory gives at least formally 

V ^ ^ J 

=0 

t 



(*) 

This integral expression is certainly a correct formula for the leading order piece of the 
wave function that made a transition after time t. However, since P^HfPjQ is of order 

one, it seems at first sight to be of order e^^~^, giving a transition probability of order 
^3^-2_ ^]^jg jg factor of larger than the expected value of order a^e~^ = e^^"^. 
Thus the integral (*) must be of order due to oscillations. We don't see any way, 
however, to evaluate (*) directly in order to get the simple formula (|9|. This is because 
the "unperturbed dynamics" given by the Born-Oppenheimer approximation is still 
a highly nontrivial Schrodinger evolution for many interacting particles. In order to 
obtain a perturbative integral expression for the leading order transitions that has less 
oscillations and is thus tractable, we replace P?o by dressed superadiabatic vacuum 
projections P|vac- Physically speaking, this is because the leading order effect of the 
coupling to the field is a dressing of the electrons of order et^ = a^. The rate of 
spontaneous emission is only of order £2^+2. However, since the dressing does not grow 
as a function of time, after times of order e'^ the spontaneous emission of order £:2^~2 
dominates the effect of the dressing. Therefore we can neglect the dressing in the final 
statement and it appears only in the proof. 

2 The main results 

In this section we only give the main theorems and explain their proofs. The more 
technical proofs of the propositions and the lemmas are provided in Sections [3] & |4] 
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respectively. In the first subsection we state a result about the time-dependent Born- 
Oppenheimer approximation and show that it remains valid, when we switch on the 
coupling to the field. In the next subsection we verify that also the superadiabatic 
subspaces survive in the coupled case. Our central results will be presented in the last 
two subsections. There we consider the transitions between different energy levels, when 



there are no free photons in the beginning. In Section 2^ we derive an expression for 
the leading order of the transition operator and in Section we provide a more explicit 
formula for the transition rate. 

2.1 Born-Oppenheimer approximation with field 

First we consider only the molecular Hamiltonian 

I 

with the usual domain -Dmoi := W^'^{R^''~^^'^) fl Hmoi- We denote the infimum of the 
spectrum of i^^oi ^ ^et A^. := Yl]=i ^^j- Furthermore, we define -D^oi := T-imoi 
and denote the maximal domain of (-f^moi)" equipped with the graph norm by D^^^ 
n e N. 

Let Ej{x) be an isolated energy band with Hei{x)Pj{x) = Ej{x)Pj{x). We will make 
use of the following version of the leading order time-dependent Born-Oppenheimer 
approximation. 

Proposition 1. Let Ej be an isolated energy band and Pj the corresponding band pro- 
jection. The operator 

H] := + (1 - Pj)HLia - Pj) 

with domain -Dmoi is self-adjoint and it holds for n = 0, 1 that 



e ^mol — e e J 



= Oie\t\+e). (10) 



Moreover, Pj G with norm bounded independently of e for all n G N. 



This is a variant of a result in |SpTe| . However, in ( 10 ) we have a slower growth of 



the bound as a function of time and better control on the domains, which is essential 



for the following. The proof given in Section |3.1| is a streamlined and improved version 
of the approach in |SpTe| . 

Now we will take the radiation field into account. Recall that 
on Dq = (Dmoi ® n (Hmoi <S) D{Hf)) and that 

was defined in (|8|. The following lemma shows, in particular, that Hq and H'^ are 
self- adjoint on Dq. 
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Lemma 1. The free Hamiltonian 



zre _2 

-"free ^ 



with domain 

is self- adjoint. The potentials Vee, V^n cind Van O'fe infinitesimally Hf^^^ hounded and HI 
and iff are Hf^^^-hounded with relative hounds independent of e. 

Hence, hy Kato-Rellich, Hq and H^ are self-adjoint on the domain Dq fore small enough. 
As before we define the "diagonal" part of Hq as 

: (P, ® 1) H'q {Pj ® 1) + ((1 - Pj) ® 1) H^ ((1 - P,) ® 1) 



-^i,field 



® 1 + 1 ® Hf. 



It is now straightforward to prove the correctness of the leading order Born-Oppenheimer 
approximation also with field, by treating the coupling to the field as a small perturba- 
tion. 

Corollary 1. For n = 0,1 and ^ < (3 < ^ it holds that 



e E 



g field 



0{e 



and 



e"'?" - e"'^^^'>firfd ) {Pj (g) 



0{e^^-'\t\+e) 



(11) 



(12) 



Proof Standard time-dependent perturbation theory yields 



because H^ — Hq = e^^ in C{Dq^^,Dq). The last statement is not completely obvious 
for n = 1, but it follows from Lemma hi and the fact that the commutator ^^^^[//^ — 
ifQ,iffj.gg] is relatively bounded by {Hf^^^"^ uniformly in e. Now (11) follows from 



= 0(e\t\+e) in £(d:;^+'.^S.oi) 

and the fact that 1 (g>e~'e^f is uniformly bounded in C{Dq) and from using the following 
technical lemma. 

Lemma 2. Let m,n G Nq with m < n. There is a constant C < oo such that if 



\B\ 



^ mol 



nJ -) — ^ /'^'^ some 5 > and a// j = 0, . . . , m, t/ien (g) 1||£(d",d™) < ^"5- 



By Proposition [T] we have Pj G and thus, again by the previous lemma, 

□ 



O 1 G This shows also (12). 
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2.2 Superadiabatic subspaces 

At first sight one might hope that the "non-adiabatic matrix elements" 

P.e-'t^>,- (13) 

from the jth to the ith electronic state are, at least for Ei < Ej, dominated by sponta- 
neous emission of photons. However, this is not the the main contribution to 



(13) comes from a velocity dependent deformation of the electronic eigenstates. More 
precisely, the range of Pj{x) is spanned by the eigenstate for a static nucleonic configu- 
ration. But the slow movement of the nuclei will deform the electronic states at order e. 



This deformation is visible in the naive "non-adiabatic matrix element" (13), but does 
not go along with emission of a photon. The projection P? = Pj + 0{e) onto the cor- 
rectly modified electronic states is the so-called superadiabatic projection associated to 
Pj. As the following proposition shows, the superadiabatic projection Pj commutes 
with -ff^oi up to errors of order uniformly on subspaces of bounded total energy. It is 
possible to go to even higher orders or exponentially small errors of the form e~e, but 
this requires some highly technical pseudo-differential calculus, cf. |MaSo2| , and is not 
needed for our analysis. 

Proposition 2. Fix an arbitrary cutoff energy E < oo. For any isolated energy band 
Ej{x) there is Eq > such that for e < Eq there are operators P? with the following 
properties: 

Pj is an orthogonal projection and for any n G No there zs > and C„ < oo such 
that for e < En the operators Pj are bounded uniformly in C^D'^^^) and 

With 1e denoting the characteristic function on the interval {~oo,E] we have further- 
more that 

\\[-^LohPj] l-B(-f^mol)||£(-^^^j^£)n ^-j < CnE^ (15) 

and 

Given two isolated bands Ei and Ej the projections Pf and Pj satisfy 

P^PnE^{HLi) = 0{E') (17) 

in ^CHmoh^mol)- 

Note that it is possible to go to even higher orders or exponentially small errors 
of the form e~e, but this requires some highly technical pseudo-differential calculus. 



cf. |MaSo2| , and is not needed for our analysis. On the other hand, the statements 



(14) and (16) without energy cutoffs and also the fact that the errors are small even in 
the norm of D^^^ for n > are essential to our analysis. Since they do not follow in 
any straightforward way from the known results, we give the proof of Proposition [2] in 
Section 13.21 
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The next natural steps would be to improve on the error in Proposition [T] for states 
with high energies cut off and to determine the asymptotic expansion of PjH^^^Pj. 
However, for our purpose it suffices to consider the leading order effective Hamiltonian 
and we will need Proposition [T] without a cutoff. Higher order results can be found e.g. 
in IMaSoillPSr;] . 

Now we investigate the lifted projectors PJ ® 1 G /^("H). Modulo some technicalities 
concerning the different graph norms, the following corollary is a simple consequence of 
Proposition |2j 

Corollary 2. Let PJ be the operator defined in Proposition^ Then for any n G No 



P= 



P,®1" 



p/® 1 



o{i), 

0{e), 
= 0{e') 



and for n > 1 



0{e) 



(18) 

(19) 
(20) 

(21) 



Proof. In view of Lemma [2|j 18) and (19) follow from the corresponding statements for 
PJ, and (21) follows from (16) and 

[HlPI®l\ = [Hl,,,P^^]®l. 

Since Hf is nonnegative and -ff^oi ® ^ ^'^d 1® Hi commute, we have that 

{1e{HIJ ® 1) 1e{HI) = lEim) ■ (22) 



Now (|20j) is a direct consequence of (15) and the simple computation 

\\[H^,,PJ®l]lEm)\\cin,Dii) = 



< 



||([i/Ll,^'] ® l)(li^(^mol) ® l)lEiH^o)\\cin,D^) 

\\i[H^^,,,P;]lEiHU ® l)lEm\\cin,D^s) 

ll[-^mob-P/]li?(-f^mol) ® M\c{DS) \\^E{Ho)\\c(n,D^)- 



□ 



In the following, we just write P? for P^ ® 1 to shorten notation. 



2.3 Spontaneous emission: transition operator 

According to Corollary [T] there are no transitions between different electronic states at 
leading order even on microscopic times of order e~^. Naive perturbation theory used 
in the proof suggests that the leading order transitions through coupling to the field 
are of order ea/^"^ on this time scale. However, we will see that the transitions through 
spontaneous emission are actually smaller, namely of order e2^~2. Since we aim at a 
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leading order expression for these transitions, we need to control the full time evolution 
up to errors which are smaller than 52/^" 2. To this end we construct superadiabatic 
subspaces corresponding to definite dressed electronic states containing so-called virtual 
but no free photons. Note that the adiabatic subspaces of Corollary [2] correspond to 
definite electronic states with arbitrary state of the field. 

For the sake of clarity we formulate our results about the dressed projector -Pj^vac 
with an additional small parameter 6 > 0. Later on we will choose a 6 that is fixed by 
e and (3. 

Proposition 3. Fix an arbitrary energy cutoff' E < 00 and let Ej be an isolated elec- 
tronic energy band. There is Eq > such that for e < Eq and any 5 > there are 
operators -Pj^vacl*^) ^ '^i'^) ^ ^{Dq) with the following properties: 
Pjvac orthogonal projections with 

II [H^ 1^+1(^0) IL(«,^o) = ^i'''^') ' (24) 

and, as a consequence, 

\\[xm,PU\\c^n,Do^ = 0{sl^6^), (25) 

for any smooth x with compact support in {—oo,E + 1). 
Moreover, for n = 0,1 

ll^- vac -P^® QoWciD^) = 0{e^^5~'^) , (26) 

where PJ is the superadiabatic projection from Proposition and Qo is the vacuum 
projection in Fock space. 

The importance of the dressed vacuum projection -Pj^vac li^s in the fact that its range 
is invariant under the full dynamics with a smaller error than PJ ® Qq. This is because 
it contains also the dressing of the electrons by virtual photons, which is crucial for 
computing the leading order transitions through spontaneous emission. The transitions 
from and into the subspace Pj^vac generated by the commutator [Pj^^^^, H^], which 
we can compute at leading order. 

Proposition 4. For < 6 < 1 and | < /3 < | 

[Pl..c,m ^E+l{Hl) = el^S^l- 1^^{H^,) + 0{el^-'\l+e5-''')) (27) 
in CiTi) with 

Tj = 15^ {Tj + T*) + 2 r ^ e{VTj + VT/) ■ eV, 

and 

T.J = -{H, + i/ei - + i5)-^Hl{P^ ® Qo) • 

For n = 0, 1 it holds 

\\'T'j\\c{D'^+\D^) = ^i^)- (28) 
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We can now compute the leading order expression for the piece of the evolution that 
makes a transition from ® Qo to Pf, which is the first of our two main results. For the 
sake of brevity we consider two non-degenerate levels Ei{x) and Ej{x) with associated 
normalized eigenfunctions (pi{x) G RanPj(x) and ^Pj{x) G RanPj(x). The results can be 
generalized to degenerate bands in a straightforward manner. 

Theorem 1 (Leading order spontaneous emission). Fix an arbitrary energy cutoff 
E < oo and let Ej and Ei be isolated electronic energy bands with Ej — Ei > and let 

I < /3 < |. Then 

lim (^£5-i/^p;^e-^^^' _ i^*e-i^^f.fleid 7-^. q-^I^U.m ds^ (pj ® Qo)1e{H') = 

(29) 

in the norm of bounded operators and uniformly on bounded time intervals. Here 

-7^- = / ^n.i'^'l'l^flt. \^.{x)){^A-)\^a*{k,X)dk (30) 

^^^'x^2,,\ ^/\k\{\k\- Ae{x) +16) 

' |A:|<Ao 

With 5 = e^-'-^-'iy^, AEix) := Ej{x) - Ei{x) and Dij{x) := Y.'e=i{'fii^)\yi'fji^))n,i- 

In Theorem |2] we will show that the integral in ( 29 ) is of order one (with a norm 
independent of the precise choice for 5(e)), although this is not obvious from the prefactor 
and the norm of Tj^i- However, this observation justifies the interpretation of the 
integral as the leading order piece of the wave function that makes a transition from 
level j to level i. 

Proof, (of Theorem [7]) In order to interchange energy cutoffs with unitary groups it 
turns out useful, that replacing energy cutoffs in terms of Hq by cutoffs in terms of H'^ 
and vice versa only adds an error of order e^^. More precisely, the graph norms induced 
by -f^o ^^^6 equivalent and we have that 

\\xm-x{m\ciH,Do) = <^i^''')- (31) 

Both claims follow from the following lemma. 

Lemma 3. Let (Hq, D{Hq)) be self-adjoint and equip D := D{Hq) with the graph norm 

II ■ II Do- IfAe C{D,n) satisfies 

\\A\\c{D,H) < 5 < 1 , 

then H = Ho + A is self-adjoint on D{H) = D and the graph-norm \\ ■ ||d^ induced by 
H on D is equivalent to \\ ■ \\oo- More precisely, for ip & D 



{i-s)\mDo < Md, < (1+5) 



Do 



Moreover, for any x ^ C'cfl^) there is a constant C < 00 (depending only on x, but not 
on H, A or 6) such that 

\\xiHo)-xiH)\\^^^^^^<C6. 
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As another useful consequence of this lemma we note that the full unitary group is 
uniformly bounded in C{Do), i.e. 

ll^~^^"1L(A.)^^ll-"'"IL(z...) = ^' (32) 

where Dffe denotes Dq equipped with the graph norm of if^. 

We now fix some x ^ C*o° with xIe = 1e and xl^+i = X and abbreviate Pj^^ := 
-^/vac- Then (using Lemma 2 from now on implicitly) 

PtPlxm ^ PtP^^®Qoxm + 0{ei^5-'^) 

^ P^P^ ®Qox{Ho) + 0{e^^{l + 5-'^))^ 0{e^H-'^) 

and, hence, 

PtPl^'^'^'P^JdH^) ^ P^P^^xiH^^'^'^'Pl + = ^(^i'^ri) . (33) 

In the following computation we make explicit the size of the error term in the line 
where it first appears and collect all of them only in the final line. We use (28), (32) 
and Pj^y G C{Dq) throughout without noting it explicitly. 
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1331 



J25l 
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|T2} 



( |3T|l,(|25t 



x{H^)PlMHn + o{el^6'^) 



(34) 



P 



e-if^' P'r 



el 
'0 



3 

1£2 



^ i£-2^-^5"^ f Pi e-^'^'^'^fi^ld 7^. e-'f^l.fleM(p. ^ Q^)x{Hl) ds 1e{H' 
Jo 



For our choice of 5 = the error term is 0(52^ 2) and we can neglect it in the 

following. 

By using that || [eV., P,] || = 0{e), P^P, = 0, T/Qo = 0, VT/Qo = and ||Vr,|| = 
(9(5~2)^ we see that 

= {ie-'H{T, + T*) + 2£5 (VT, + VT/) ■ ^V.) (P,- ® Qo) 
= ie-'26PiT,{Pj ® Qo) + 2e^PiVTj{P, ® Qo) ■ ^V, + 0{eh 



3 ^ 3 , 
2 
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Thus the corresponding replacement in the integrand contributes an error of order 
C)(e2^+2(5~2 |t|). Since V^-Pj and V xPj ® Qo are bounded independently of e and, since 

II (iff + H,,{x) - E,{x) + i5)-^\\ < , 

we have 

£^P,VT,(P,®Qo) 
= -£^P,V {{Hi + H,i{x) - Ej{x) + i5)-'Hl{P, ® Qo)) {Pj ® Qo) 
= -£5V {Pi{x){Hi + ifei(x) - Ej{x) + i5)-^Hl) {Pj ®Q^) + [e^H^^^ 

= -e'^PiV {{Hf + Ei{x) - Ej{x) + iS)~'Hl) (P, ^ Qo) + O (^^<^~') 

= e"^Pi{VEi{x)-VEj{x)){Hi + Ei{x)-E,{x)+i6)^^H'^{Pj^Qo) + O (eh- 
in the norm of bounded operators. Recall that 



i=i 



3 

and thus on the one particle sector of Fock space, abbreviating p{k) := (27r)~2 l[Q^^y](fc), 
PiHl{P,®Qo){x,k,X) = 

= -4v^E -^\v^{x)){^.{x),e-"'■y^e,{k) ■ V,,^,(x))^,(^,(x)| 
i=i v^l^l 



t V2|A;| 



2v/^(E,(x) - E,{x)) -^j=ex{k) ■ J2{v^ix),ymix))nJv>^ix)){VJ{x)\ + 0{e 



m 



1=1 



27r {E,{x) - E^{x)) ^e,{k) ■ D,,{x)\v,{x)){v,{x)\ + 

\k\ 



Collecting the previous observations, we showed (29) with a transition operator given 

by 



i6 



27r£2 



A=l,2 



exik)-Dij{x) Ae{ 



X) 



\k\<Ao 



\k\{\k\-AE{x) + i5) 



X 



2ieV,As(x) 



'-i£:V.x) X 



6^\k\{\k\-AE{x)+i6) 
\<f^{x)){y^j{x)\ ® a*{k) dk + O (^£^^5 + eh-^ 

--: h + t2 +o(£^-^ + £^ri) . 



(35) 
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The error term is o(l) and we can neglect it. While the norm of t2 is of order (p-)^ and 
thus slightly smaller than that of ti, which is of order (-)^, both grow as e — )■ 0. In order 



to show that the contribution of ti to the integral in (29) is of order one and that of ^2 
is strictly smaller, one has to perform the time integration and use cancellations due to 
oscillations. In order not to duplicate the corresponding arguments, we skip the proof 
that the contribution of ^2 is negligible at this point and comment on it instead after 
the proof of Theorem [2j At that point we will have introduced the necessary machinery 
in order to explain the argument. □ 

2.4 Spontaneous emission: rate of decay 

In order to obtain an explicit formula for the leading order rate of spontaneous decay 
from RanPJ to RanPf we evaluate the norm of the leading order wave function in 
Theorem [U 

Theorem 2 (Probability for spontaneous decay). Under the same hypotheses as 
in Theorem\^for 

it holds that 

^ra(^e^-'^me--^l^'n^^ = (36) 

uniformly on compact time intervals and uniformly in ip. 

Before we come to the proof, we collect some remarks on the result: 



1. Note that the subtracted term in (|36|) is of order 1 and therefore the same is true 

;i we 



for e^~^f^ ||P/e~'e^^ Now recall that £^ is equal to the coupling constant a 
by choice of /3. So ||P/e~'e'^^ is proportional to and grows linearly in time. 
Observing that in our units time is scaled with a^, we see that Theorem [2] is the 
generalization to molecules of the physics textbook result that for atoms the decay 
rate is | AiPA|. 

2. For /3 < 1 we may replace Pf and PJ by Pj and Pj respectively in the theorem 
because — Pk = 0{e) for k = i,j and e^^~^ > e. 



3. Practically the result (36) means that the decay rate can be computed at leading 
order within the Born-Oppenheimer approximation: while i\} G "Hmoi still contains 
the electronic degrees of freedom, Pji\} is of the form (f){x)(pj{x,y) and therefore 



(e-^t^lP,V^)(x,y) = (e-t^V)(a:)¥^,(x,l/) 



where the effective Born-Oppenheimer Hamiltonian is just 
hi, = (-iV. + Aj{x)f + E,{x) + 0{e^ 
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where Aj{x) := i{ipj{x), V(pj{x))y,^^ the connection coefficient of the Berry-connec- 
tion V^. := Pj'VxPj- Hence the decay rate can be written as 



||||A,|Af e-'f-IP,^ 



and one only needs to solve an effective Schrodinger equation for the nuclei. 
Proof, {of Theorem^ According to Theorem [l] we have that 

lim (^^-^/^lli^^e-'i^^ - W) = 

with 



0,,(t):= ];;dre^?^^f7-^,e-'i^):fv^ 



(37) 



and we need to show that lime_j.o (^^ijit) — | || |Z^jj|A^^^ e dsj = 0. First 

note that because of the explicit form (30) of Tj^i this state lives only the one-particle 
sector of Fock space. Writing k = u\k\, the only dependence on the angular variable u 
appears in the polarization vectors ex{uj). Using that for any symmetric 3 x 3 -matrix 
A it holds that 

Ea=i,2 Is^ do; {exiu), Ae^iu)) = f trA , 



one can perform the angular integration in (37) and obtains a factor However, to 



not overburden notation, we will make this exp 
as 



icit only later on. First we rewrite (37) 



' l^ds r^\s' e'(^+^)^^.*7-%,e-'^'^?.f7-_,e-^(^-^)^^>^vi/ 







a{s) 



(38) 



=:/(s,s') 



where s = {r + r')/{2e), s' = {r — r')/e, a{s) := min{2s, 2(- — s)}. Let 



Hj := P,i-e^Ax)Pj + Pi{-e^Ax)P, + 



and 



Hi := P,i-e^Ax)Pj + Pi{-6^Ax)P, + 

which satisfy 

PjHj = HjPj = PjH]Pj , PiHi = HiPi = PiH^Pi and Hj - Hi = A . 
Then with Tj^i = PiTj^iPj and = if) ®Vt we have that 



^, e^(^+^)^IP,-7^%,Pie-^^'^aPi7^^iP,e-'(^-1^)^l^ 
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Making the angular integration explicit, we can thus replace I{s, s') in (38) by 

' 3 i2ny e Jo V' R - Aix) - i6 

^ ^ ' R - A(x) +i6 ^ 

= :U{s,s') ^ ' 

^ i?-A(s + ^)-i(5 ^ ^ R-A{s-'^)+i6 



Sire 







Here and in the following we denote for any operator O the Heisenberg operator e'^^^ O e"''^^^ 
by 0{s) and for better readability we abbreviate R := \k\ and D := Dij. Moreover, we 
will still write out 6 in the expressions and in most remainder estimates, but keep in 
mind that in the end we put 6 = e2~^^~i^/^. 

Next we show by a stationary phase argument that / is small for large s'. Integration 
by parts in 



shows that 

/ = 0{6-^6-^T-^ + S^e-^T-^) for > r and r > . (39) 

Instead of giving the detailed computation we just mention that the boundary terms 
contain the operators 

and 



-A(s±|)±i5 Ao-A(s±|)±i5 

which are all uniformly bounded, since < A(a;) < Aq uniformly in x. So the first 
boundary term, which is of order 6'^e~^T~^, is indeed the worst. 

Next, for |s'| < r we expand the operators around s' = 0. Clearly 

A(.±^) = A(.) ± i A(.±^)]d." = A{s) + 0{eT) (40) 

in C{Dq, %) because the gradient of A is bounded independently of e. The same is true 
for A replaced by D because y?, and ipj as well as their derivatives with respect to x 
decay exponentially in y (the proof in |WaTej is easily adapted to unbounded potentials 
whose derivatives with respect to x are bounded). By using the so-called Strang splitting 
(see |JaLu] ) we see that 

= Ss+i)H^^4H^^s'Ai.)^-4H^^-iis-^)H^ ^ 0(r3||[A,[iJ|,A]]||) 
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in C{Do, %). Plugging these expansions of A, D, and U into /, which is allowed because 
^ & Dq due to the energy cutoff and all operators involved are in C{Do) with a norm 
bounded independently of e, 6, s and s', we find that for |s'| < r 



^ '"in - A(s + - W R- A(s - + W 



2^ " " V" 2 
V 



=:J{s,s',R,5) 

+ 0{t + er^) . (41) 
While this is not small, note that the contribution of the error term to Bjj is of order 



which is indeed small for our choice of r = e~5+('^~i)/^o. 



2_4^ 



When expanding the denominators in (41), we have to be more careful, since e.g. 
1 1 



i?-A(s + f)-i5 i?-A(s)-i5 



(42) 



^ -(A(. + ^)-A(.)) ^ 



i?-A(s + |)-i(5 ' ' " ' R-A{s)-i6 



is only 0{^) when naively estimating the norm, which yields a term of order 0{^) 
to Qij. This will be large for our choice of 6{6) and t{6) and thus we need a better 
estimate. For this we have to evaluate the integral explicitly. In order to prepare for the 
residue calculus, we first show that we can extend the /^-integration to all of M with a 
negligible error. Adding the integral to +oo yields 



25^ 

3ne 



[ eds' Mm dR Re-'''^J{s, s' , R, 5) 

J-T Jao 

lim f eds' I dR (i£e-'''^yj{s,s',R,S) 

^"^"""^ J —T J An 



lim f eds' r dRe~'''^-^J(s,s',R,6) 

+ l^e lim [ dR (e-'^^J(s, r, R, S) - e'^^Jis, -r, R, S)) 

= 0{6't + 6') ,° (43) 

since < A(x) < Aq uniformly in x implies that J{s, s' , R, 6) and ^J{s,s',R,5) are 
uniformly bounded by in C{Dq,1-L) for i? > Aq and analogously for i? < 0. Thus we 
can integrate from — oo to oo in i? while adding an error of order to Ojj. 
Now we want to replace J(s, s', i?, 5) in (41 ) by 



J{s, s' , R,6) := — -— —e 



R-A{s)-i6 R-A{s)+i6 
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One of the two terms appearing in the difference is 

1 

R-A+-i6 ~ R-Ao - i(5 y i? - A_ + i(5 



hm / ' dRRe-'^ - — ^ ) f{s)——— , (44) 



where we abbreviate 



A±:=A(s±^), Ao:=A(s) and f{s):=D*{s)D{s)A{s)\'''^^'K 

To show that this term (and analogously the other one) gives only a negligible contribu- 
tion to Qij, we use the residue calculus. For s' < we need to close the contour in the 
upper complex plane. Writing the spectral representation of the self-adjoint operators 
A-t resp. Ao with spectrum contained in (0, Aq) as 

/■Ao 

A+n=: / XdPt^^ 







^±,0 =• / ^U/^A ' 

the residue theorem yields for s' < 



p 

lim / dRRe-'"'"^ ( ^ - f ) /(s 



1 1 \ 1 



3'^^ p'To^ J \R-A+-i6 R- Ao-i6 J R- A_ + i6 

-p 

Ao P 

hm / / dRRe-'^'^ ( -^^^ - „ ^f" /(.) ^ 



p':;'^J J _ A - 15 R - X-i6 J R- A_ + i6 

-p 

Ao 

l£/Ae-^'(^-'^)(dP,+ -dP°)/(.)3^-^l^ + 0(f) = (*). 







According to (42) we can replace A_ by Aq in the resolvent at the price of a term of 
order — ■ p- = t. Then one can commute the resolvent with f{s) and afterwards replace 
by the same reasoning Aq by A4. for the first summand. Now one can integrate the 
spectral measures explicitly again and obtains 

Aq 

U) = ^ [ Ae-'^'(^+'^) ( ^ ^ 1 f(s) + 0('- + t) 

7 VA-A+ + 2i5 A-Ao + 2i5;^^^ Ve^J 



/„-is'A+A „-is'A() A \ 



0(tH+'-^+t), (45) 



where we used that by exactly the same reasoning as in (40) we have 

e-'^'^+A+ - e-'"'^°Ao = 0{eT^) 

for < r. The additional factor of r comes from the fact that derivatives of e^'^'^^^^ are 
of order \s'\ < r. After integration over s and s' this adds an error of order er^d+rd'^+eT'^ 
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to Qij. To estimate (44) for s' > 0, one closes the contour in the lower complex plane 
and proceeds along the same lines as above. 

Let ar{s) := min{a(s),r}, then collecting once more all the estimates we obtain 



_2 



/ ds / ds'J(s,s') / ds / ds' lim / dRRe-''^ J{s,s',R,S) 

Jo J -I Stt Jo J-a,(s) p^°°J-p 



P9l 



^ +8^ dsj^ ds' O [t^6 + f + 
= O {{5T)-'e-' + d'e-'r-' + 5't + + er^ + e'r^ + . 

Recall that 6 = and r = ^-l+C/^-D/io and thus 6t = £-(/5-i)/io ^ i fo^ 

I < ^ |- Then for / big enough, the error is o(l) for all /3 with | < /9 < |. 

Finally, we compute the main term using again the residue calculus. We close the 
integral depending on the sign of s' and get 



2pe 



25^ 
37r 



25e 
3 



2fe 
3 



/ ds ds' lim / dRRe-''^Jis,s',R,5) = 

Jo J-ar{s) P^°°J-p 

t'ds r'ds' hm rdi?i?e-'-^^:^M^e-'-W^M^ 
/'ds P^'^ ds'e-l^'l'^|D|2(s)A3(s) + 0{5^t) 

Jo J-ar{s) 

P - POO 

I ds / ds' e-l^'l-^ |Dp(s) A3(s) + C((5V + e""-^ + er) 

Jo J-oo 

t 

f ds\DWs)A%s) + 0{5^T + e-^^ + eT). 
Jo 



So we end up with 

fa(s) 



0^^.(t) = £3/3-1/^^^2 ds'/(s,s')^ 
\ Jo J -a(s) I 

= e^^ r^'|(^, |/^r(s)A=^(s)^) ds + o{e^^-') 
Jo 

= e'^-' r||||A,|Af e-^?^lP,V^||'ds + oie'f'-'] 
Jo 



□ 
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We still need to show that the contribution of t2 from (35) to the transitions is 



negligible at leading order. More precisely, we need to show that 



Qi^it) := [ dr [ dr'e'?^I.ft;e-'?^-ve'i^^'.ft2e-'^^^>f^' 

\ Jo Jo I 

v]/, f ds r'^dsV(^+1^)^^.^t;e-^^'^^,^t2e-'(^-^)^^^f^\, 

Jo J-a(s) ' 



is oi^e^^ With the same type of arguments as in the previous proof one can now show 
that the main contribution to Qij{t) comes from the integral 

"Jo \ J-r P-*^J-f (JS - A(s) - 1^)2 (R-A(s) + iSy 

which is easily seen to be O (^e^ (^-^ + j after performing the R integration. Here 
^e{s) := e'^'^^-iedcje''''^^'^ E Dq. This concludes to proof of Theorem [l] 

3 Proofs of the main propositions 

Before giving the details of the proofs let us shortly comment on the relation and dif- 
ferences between the Propositions [T} [2] and 3. In some sense they are all "adiabatic 
theorems", however, of slightly different spirit. In Proposition [T] we adapt and simplify 
arguments from |SpTe| , which in turn were motivated by Kato's proof of the adiabatic 
theorem of quantum mechanics for Hamiltonians slowly depending on time. The basic 
idea is to show that the transitions between adiabatic subspaces are small even for long 
times by explicitly evaluating an oscillatory integral. In Proposition [2] we use the idea 
of superadiabatic perturbation theory: the adiabatic subspaces are replaced by slightly 
tilted superadiabatic subspaces. The coupling between the superadiabatic subspaces is 
so small that the transitions between them can be estimated even for long times by a 
crude norm-estimate of the integrand. The technical reason that forces us to include the 
weaker statement of Proposition [T] is that it can be easily proven without energy cutoffs. 
This is crucial when replacing the full time-evolution by its adiabatic approximation in 
the computation (34) in the proof of Theorem [l} 

An essential input for adiabatic decoupling and thus for all proofs in this section is 
the fact that the smoothness of He\{x) and the gap assumption imply the smoothness 
of the map Pj : M?^ CiUd), x ^ Pj{x). 

Lemma 4. Let Ej be an isolated energy band and Pj the corresponding band projection. 
Then Pj e (M^', £(Hci)) and Ej e C{^(M3'). Moreover, for any a e N^^ and any 
n G No one has d^Pj{x) G C{HehD{Hl)). 

Apart from Pj{x) there will appear numerous operator- valued multiplication oper- 
ators of this type and in addition also differential operators 5^ with operator-valued 
coefficients. The following lemma will turn out useful when working with these kind of 
operators. 
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Lemma 5. Let '■ M'^' — )■ C{'Hc\) be hounded, smooth and with hounded derivatives, 
i.e. Aa G {E?\ C{T-Le\)) . Then Aa = Aa{x) dx defines a hounded operator on 
^moi = L'^{^^'']'Hei) and we call the differential operator 

n 

A^=Y. 

|a|=0 

an admissihle operator of order n. 

(i) An admissihle operator of order n is hounded in £(Z)JJJ'qj, "Hmoi) for m = [n/2] 
uniformly in e > 0. 

(ii) // all coefficients Aa of an admissihle operator A^ of order n have the property 
that [{H^^^)'' , Aa] is an admissihle operator of order 2k — 1, then is uniformly 
hounded m C{D^J^, D^^^) for m = \n/2] . 

As a first simple application we note the following corollary. 

Corollary 3. d^Pj is uniformly hounded in C{D^^y) fof^ G No and (3 G Nq'. 

Proof. According to Lemma |4| 9f Pj is an admissible operator of order for any (3 G 
Nq'. Statement (ii) of Lemma [s] implies the claim of the corollary once we show that 
[(-^moi)") -^i] admissible operators of order 2n — 1 for any n G N. This in turn 
follows from direct computation and the fact that according to Lemma |4] we have 
d-P,ix)eC{n,i,DiH:,)). □ 



3.1 Proof of Proposition [T] 

Since ed^^ has norm one in C{D^^l , D'^^^) , Corollary ^ implies that the commutator 



-e\A,Pj 



2eW^Pj ■ eV, 



0{e) 



(46) 



is of order e in CiD"^^], D"^^^) for all n G Nq. Set P/(x) := 1 - Pj{x). Since 



TTS TT£ 

-"mol -"j 



;i-2P,)[i7^„i,p,], 



when 



the self-adjointness of (f/J, -Dmoi) for e small enough follows from Lemma [3j 

We notice that ||e~^£^m°i = 1 and \\c{d{h^)) = 1 for all t G 

D{Hj) is equipped with the graph norm. Then the equivalence of the graph norms due 

to Lemma sjimplies that ||e~'^^J ||£(Dmoi) bounded independently of e. 
Now set Rj{x) := P^{x){H^x{x) -^Ej{x))-'^P^{x) and 



K,{x) 



R,{x)Hl,,P,{x) + P^{x)HlM^) 
R,{x)[Hl^,,P,{x)]P,{x) - P,{x)[Hl,,,P,{x)]R,{x) 



Due to (46) we have that [H^^y, Pj]Kj, [Ej,Kj], and [—e'^Ax,Kj] are of order in 
C^moi)- Thus it holds that 

[i/Li, K,] = [//el, K,] + Oie') = P^Hl,,P, + P,Hl,,Pf- + 0{e^) (47) 
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and therefore 



, .e ^ ^ K^e 'e^mol 

as 



= ^e-'^^^l((i/Li + (2P, - 1)K,,,P,])K, - K.nQe- 
in -Dmoi)- Hence the difference in the unitary groups is 



e mol 



e 



^0 



'0 

n+l Y)n 



Since i^j is of order e in , -DJJ^oi), we obtain (10). 



3.2 Proof of Proposition [2] 

This construction has been done in different places using different techniques. For the 
most general treatment of the Born-Oppenheimer approximation allowing even nuclei 
that are point charges we refer to the recent work of Martinez and Sordoni |MaSo2 



based on a twisted pseudo-differential calculus. Since the precise statements we need 
for treating the coupling to the field do not follow from their results, we give a more 
elementary proof for the case of smeared nuclei here. It is partly an adaption of the 
arguments used in |WaTej in a different context. 

For better readability we now drop the index j and write Pq •= Pj ^^nd := Ej. 
To have some margin to play with we use first the characteristic function 1e+i on 
{—oo,E + 1] and recall that with e denoting the infimum of the spectrum of -f^moi 
have that 1e+i{H^J = l[e,E+i]{H^J. 

Starting from the orthogonal projection Pq we want to construct a self-adjoint oper- 
ator P*^ G CiJ-Lraoi) with 

P^P^ = P^ and P^] = Oie') . 

The first statement just means that P^ is a projection. The basic idea for constructing 
P^ is to determine first the coefficients in an asymptotic expansion 

P" = Po + ePi + e^P2 + 0{e^) , (48) 
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where we recall that according to (46) the commutator [H^^^,Pq] with the choice Pq = 
Pj is of order e as an operator in DJ^^j). As shown in many instances, the 

requirements that P^'^^ := Pq + ePi + e^P2 satisfies 

p(2)p(2) _ p(2) ^ 0^^3) [P(2), i/Ll] = 0{e') 

fix P*^^) uniquely modulo terms of order e^. We will not repeat the construction here, 
but only give the result: Let 

[-Po] := \ [H^oh -Po] , 

then [Pq] is, according to (46), a uniformly bounded operator in C{D^^^ -^moi)- We put 

:= Po[Po]P , 

with the reduced resolvent R{x) = Po{x)^{Hci{x) — E^:{x))~^Po{x)^. Since d^Hei{x) 
and thus also d^R{x) are bounded operator on Hei for any a G N'^', 5*1 is an admissible 
operator of order one in the sense of Lemma |5] This is where smearing out the nuclear 
charge distribution is essential. By the same reasoning as in (47) this choice makes 

[Po] + [H,y,S, + Sl]=0{e). 

Now let 

Pi := ^1 + ^1 and P^^^ := Pq + ePi , 

then 

p{i)p(i) _ p(i) = s{PoPi + PiPo- Pi) +e^PiPi 

= eiSi + SI - Pi) + e^SiSl + SlSi) = e^SiS*, + S^Si) 

and 

[i/^i, P(i)] = e[Po] + e[H,,, Si + S*,] + e^[-eA, Pi] = 0{e'') 

in C{D^l, ^moi)- Now we simply iterate this construction: first we modify P^^^ in order 
to make it a projection to higher order by putting 

p(i) pW + e^SlSi - SiSl) . 

This gives 

p(i)p(i) _ p(i) ^ p(i)p(i) _ p(i) _ 2e^SiSl - e\SlSi - SiS*,) + 0{e^) = 0{e^) . 
Then we put 

and 

which makes 

[P^'^] + [H,i,S2 + S;] = Oie). 
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Defining 

Jr^2 — ^2 ~r 02 -r OiOi — OiO| 

we find that 



and still 

p(2)p(2) _ p(2) ^ p(l)p(l)+^2^^^^^*)_(p(l)^^2^^^^^*))^(>^(^3) ^. ^3pe_ 

Note that i?^ is an admissible operator of order three and i?| is an admissible opera- 
tor of order four. The following lemma shows that all the operators appearing in the 
construction can be bounded by appropriate powers of -f^^oi- 

Lemma 6. The operators Pq, Pi, P2, R\ and R2 are admissible operators of order 0, 
1, 2, 3, and 4 respectively. Their coefficients have commutators with (-f^^oi)'^ ^^^^ ^''"^ 
admissible operators of order 2k — 1 for any /c G N. Thus they are uniformly bounded 
operators from to D'^^y /^^ ^''^V n E'Rq and m = 0, 1, 2, 3, and 4 respectively. 



In order to make sense of (48) as a bounded operator in £(7^moi) and to get uniform 
bounds on [-ff^oi; t\ms need to cut off large energies. To do so we fix < 00 and 

choose xe € Co°(IR, [0, 1]) such that XE\[e-i,E+i] = 1 and suppxs C {e — 2,E + 2). Then 
we define 

P' ■=ePi+e^P2 

and 

P^^ := Po + eP, + e'P2 - (1 - Xe{HI,,)) P' (1 - Xe{HI,,)) , 

i.e. we cut off the corrections to Pq at high energies. To see that P^^ is indeed a bounded 
operator in C{D'^^^ for any G No? note that 

P^^^ = Po + P'xe{HIJ + XE{Ht^.;}P'{l - Xe{HIJ) 

and that Pq and P^xe{H^^{) are bounded independently of e in C{D^^i) by Lemma [g] 
and the fact that Xe{H^^^) G /^(Hmoi, D^^^) with norm bounded independently of e. In 
particular we have also 

\\PL-Po\\^^^r.j = Oie). 

We first proof that the operator P^^ has all the properties claimed in the proposition 
modulo the fact that it is not a projection. In a second step we turn it into a projection 
without loosing the desired properties. 
Now by Lemma |6] it follows that 

as a bounded operator from D^^l to D^^^^. With Xe{H^oi) ls+i(-f^moi) = '^e+i{H^oi) 
this implies 
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as an operator from Timoi to i^^oi- Note that, by taking adjoints, this imphes that 
and with ||lE+i(i^^oi)lk(Hw,^^„i) = ^(1) also 

For later use we also show that this implies the smallness of the commutator of P^^ 
with a smooth energy cutoff x supported in (e — |, + |) 

ll[x(a-..).-p«]IL,«„.„z>,.,, = o(=')- m 

Since the argument will be used several times in the remainder of the paper, we formulate 
it as a lemma. 

Lemma 7. Let I C M. be a compact interval, I d I another interval with different 
endpoints and x G C'^(M) with suppx C I. Then for any n & Nq there exists C < oo 
depending only on n and x with the following property: Let {H, D{H)) he self-adjoint 
and A e C{H) he hounded and self-adjoint. Then 

\\[H,A]li{H)\\c(H,DiH-))<5 

implies that 

ll[x(^),^]ll£(«,i.(H"))<^^^- 

Now we need to turn the "almost projection" P^^ into a true projection. Since we 
will use this trick as well several times, we formulate it again as a lemma. 

Lemma 8. There are constants Cn < oo, n e N, such that the following holds: 

Let {H, D{H)) he a self-adjoint operator and let := D{H"') he equipped with the norm 

n 

:=5^||i/V||- 

i=0 

For some N E N let Q be an operator that is hounded in C{D'^) for all < n < N and 
self-adjoint in C{l-L) with the following properties: 

WQQ - Q\\c(D^) < S (50) 

for all < n < N and some S < ^ and 

\\[H,Q]\\ciDn,Dn-i)<5n (51) 

for all 1 < n < N and some dn < 

Then there is an orthogonal projection Q G €,{%) with \\Q\\c{D") < ^^'^^ that satisfies 

\\Q-Q\\c[H)<5 and \\Q-Q\\c{Dr^)<Cn5 (52) 
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and 

m,Q]\\ 

for all n < N . 

Moreover, there is a constant Ce depending only on E E M. such that we have the 
following implications: 

\\[H,Q]lE+iiH)\\c(^n,D-)<(3i (53) 

for all n < N implies 



for all n < N, and 
for all n < N implies 

for all n < N . 



\\[H,Q] lE^.^{H)\\a^n.D-)<CECnPi 

\\{QQ-Q)lE+\{H)\\m,Dr^)<P2 (54) 

m - Q)lE+dH)\\m,Dr^) <CECnP2 (55) 



We can now apply Lemma ISj to the almost projection P^^ almost commuting with 
if^Qi, where now 5 and 5„ are oi order e and (5i of order e^. To be able to use also the 
last implication of Lemma [s] with of order e^, we still need to show (54). To this end 
observe that we have by construction that 

||(p(2)p(2) _ p(2)) 1^^, = 0{e') . 

Hence for x ^ C'^(M) with xl^+f — suppx C (e — |, + |) we have 

= 0{e'). 

Thus we can use Lemma|8]to turn P^^ into an orthogonal projection with the desired 
properties. 

Next we show that 

pj^)pf ) = 0{e^) for t^j. 

To enhance readability we denote Pj^^ = P^^^ and Pj^^'' = Q'-^-' etc., i.e. we distinguish 
the different electronic levels by the letters P and Q instead of the indices j and i. Then 
with Qo^o = 

gWpd) = (Qo + eQ,){Po + ePi) = e{QiPo + QoPi) + 0{e^) . 
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Denoting the S'l-operator associated to Q^^^ by Ri we find 

QiPq + QqPi = RiPq + QqSI 

= Qo[Qo]RiEQ)Po-QoR{Ep)[Po]Po 

= {Ep-Eq)-\Qo[Qo]Po + Qo[Po]Po)=0. 

This imphes that 

P^gn^^.(i/Li) = P^Q^xiHWi-E+iiH^moi) 

= 0{e') 



in /^("Hmoi, 



mol/ 



3.3 Proof of Propositions |3| & g] 

We first recall the perturbative form of from ([S]): 

=: + e^^Hl + ei^+^i^l . 
The operators Hf and i/| satisfy 

\\Hl\\c{D{,H^,).M) < 

with constants Ci independent of e. Hence Lemma [3] yields that for e small enough 
is self-adjoint on D[Hq) and the graph norms induced by Hq and are uniformly 
equivalent. 

We write P^^^ := -Pj^vac ^'^d Pvac := -Pj ® Qo, where Qq is the projection onto the 
vacuum state in As before we first construct an almost projection P4,c with the 
desired properties and then apply Lemma [8j Since the first correction to Hq is of order 
£2^^ it is natural to make the ansatz 

P4, := P^®Qo + el^Psp , 

where we assume that P^ is constructed as in Proposition [2] but with energy cut off at 
E + 1. Computing the commutator with if^, we find that 



(l56l r^. ^ ^ . 

2 



J57b-d59l 



[P^ ®Qo + ei^Pa^, H^, + e^^Hl] 1e+i{H^o) + 0{elf'+') 
[P^ ® Qo, ^0^] lE+im + [Pp. HI] li^+i(^o) 
+ ([P^®Q,,Hl] + [P.^p,Hl]) lE+im) + Oiel^'^') 
e'^' {[P...,Hl] + [Pi„H^,]) 1e+,{H^o) + 0{ei^^') 
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in In this computation we made the following assumptions, which are clear on a 

formal level but need to be proved later on: 

W(i7o^) = 0(l) (56) 
e'^[P.,,Hl]lE+,{Hl) = 0{e'^'^') (57) 

[P^ ® Qo, H'o] li?+i(^o) = 0{el^+') (58) 
[(P^ - Pj) ® Qo, HI]1e+i{HI) = Oie) , (59) 
all in £{1-1). Whether (56) and (57) are satisfied depends on Ps^, which we now construct 
by the requirement that the commutator is of order e^^S'i, i.e. that 

Dropping the energy cutoff for a moment this translates to 

[P3^,i/0'] = -[P.ac,/^l^] + O(55). (60) 

To solve this equation for Pa^ we cannot proceed as in adiabatic theory with spectral 
gap, since the reduced resolvent (ifei + H{ — Ej)^^{l — Pvac) is not bounded without 
spectral gap. Therefore we proceed as in |Teui| and shift the resolvent into the complex 
plane by a small amount 6, 

Pip := -{H, + H,,-E^ + i5)-^HlP,.,,~P,^Hl{H, + H,,-Ej-i5)-' 

Note that Pll,^ is exactly the first order correction one would obtain by formally applying 

standard perturbation theory to the electronic eigenprojection Pj{x) ® Qo- With this 
definition we find that 

[P|^,iff + iJel] (61) 

= HIP,^ - P^^M + {Ej-i5){H, + H,,-Ej+i5)-'HlP,^ 

- P..,Ml{Hi + i/ei - E, - i5y\E, + i5) 

- (Hf + H,, - E, + i5)-^HlP,^,E, + E,P,^,Hl{H, + H,, - E, - i5)-^ 
= - [P,^,Hl] + i5{Ts + TI) . 

We will show that 

i5{Ts + T;) = 0{5"') and [T,, eV.] = O (^) , 

which indeed gives us (60) for (5 > e^. Note for the following that P|„ is again a fibered 
operator. 



*mol 



This is important, since we will need to commute P| through H^^^ and thus to compute 



derivatives of P3 (x) with respect to x. 
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Lemma 9. For 6 > small enough we have that M^' — )■ x t— )■ Ts{x) is smooth 

and there is a constant C < oo not depending on 5 or e such that 



< 



C_ 

71 



\\{H, + H,,)Ts\\ < 



C_ 

71 



(62) 
(63) 



and for \a\ E Nq' 



Wdtm < c 



H + 5 



and \\{Hi + H,{)d'^Ts\\<C 



|a| + i 



(64) 



Proof. Let {ipi{x), . . . ,ips{x)) be an orthonormal basis of RanPj(a;) and write \1/ G 
RanPj ® Qo as 

s 



Then 



m=l 



s r 2 



//[Pvac* = iJ2Me'y^) ■ V..^ = ij] 5^X^^„(a:)^p(A:)e-''=-^'V,,V^^(a:,y) 



i=l 



m=l i=l A=l 



Since the sum is finite, it suffices to estimate the resolvent acting on each summand. 
We spht 

1 = l[e,oo)(-f^el) = 1 [£,£.] (-f^el) + 1 (E, ,00) (-f^cl) =: P< + P> 

and observe that on the range of P> the resolvent is indeed uniformly bounded also for 
5 = because of the gap condition. So it remains to look at the resolvent acting on the 
range of 

£=1 

Using Hp]Pf = EfPf we get that 



+ H,,{xmk\ + H,,{x) - E,{x) + i5)-^-^P<{x)e-'''^-y^Vy^^m{x, y) 
j 

J2{\k\ + E,{x)){Ee{x) - E,{x) + \k\ + 16) 



-1 m 

c 

< _ 

^Jo {\k\-{E,{x)~E,{x))f + 5^ - 5 



P,{x)e-'''^-y^Vy^^m{x,y) 



This shows (62) and (63). 
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To get the bounds for the derivatives first observe that whenever a derivative hits a 
resolvent, we get 

d^^{\k\+H,,{x)-Ej{x) + i5)-^ 

+ H,,{x) - Ej{x) + i5y' d.^Ejix) - Helix)) {\k\ + H,i{x) - Ej{x) + i^)"' 



where dx^Ej^x) — Hci{x)) is uniformly bounded. By Lemma |4] derivatives of Pj map 
into the domain of ifei and thus into the domain of iff. Hence, whenever at least 
one derivative hits Pj there will be at most |a| resolvents left and such a term can be 
estimated by 5~I"L When all the derivatives hit the resolvent, the worst term has \a\ + 1 
resolvents, which can be estimated by times the norm of T^. □ 

Corollary 4. Let < 6 < 1. With 

Py ■= Ts + Ti 

we have that for n = 0, 1 



\pi 



p\\c{D{{H§y-)) = 0{S 2), 



and 



Moreover, 



where 



I pe p£ 

I vac j 



l-^vac-^vac ~ Pvac\\ciD{iH^)")) 



0{e 



PL.H, 



) -'-'0 



+ [Po®Qo,m] = 5-^T + R, 



(65) 

(66) 
(67) 

(68) 



r = i6-^ (Ts + r;) + 2 r 5£v.(T5 + r;) ■ ev. = o(i) 

zn £(D((ifg)"+i),L)((/7g)")) and 



\R\ 



Proof. We will use 6 > without noting it explicitly. It follows directly from (62) that 
IIT^II = ||T;|| = C((5-5), which yields Q forn = 0. Forn = 1 notice that (Q, 
and Ts{Hei + Hi) = TsEj imply that 

[H',, Ts] = [-6^ Ax + H,i + Hi, Ts] = [s^Ax, Ts] - H^P,^ + iSTs 



is C(l) in C{D{H^),n) and thus by Lemma [l| also Ts = 0{6-^) in C{D{H'q)). In 
a similar way we find that also T; is 0(5-2) in C{D{H^)). The estimate m§ follows 



immediately from (65). For (67) we note that 

Pva.cTs = = 



rri* p 
S vac 



(69) 
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because [Pvac, -f^ei + H{] 
have 



and Hf creates a photon when apphed to Pvac- Then we 



2 2 



pe 



Qo + Po®QoPL)+0{e'^5-') 



J69l 



For (68) note that according to (61) 



With 



and IIP] 



PL 



Hp] + Hf 



[Po®Qo,Hl]=i5{T + T*). 



cm 



\e\A,Pf 



ii3)\\c{n) 



■-2e{V^Pl)-eV. + e'{A,P! 

2^ 2 

(9(5^(5^5) we directly obtain (|68|) for n = 0. For 



n = 1 it suffices to show || [T, -f^ol = It holds 



As shown in the proof of (65) the ffist term is of order 62. Analogously, it follows that 
the second term is of order e^"^. Hence, both are 0{1) because of < 6 < 1. Finally 

\C{D{H^)) = \\e^{A^PL)\\c{D{H^)) 



\R\\an(m)) = ll^"^ (A^PI o)\\c(d(hf.)) = 0{e^8 2+^2^ 2) follows from Lemma 



and (l64j). 

Lemma 10. It holds that 



12 



□ 



\[H'.Hj\\ciDiH^,)M) = 0{e + el^5'^) 



and 



[P. 



in OT-L). As a consequence, 



xm,p. 



e 

vac 



c{n) 



(70) 

(71) 

(72) 



for any smooth x with compact support in {—00, E + ^] 



Proof. We ffist show ( 71 ) . Due to (|68| 



we only need to check (56)-(59). The estimates 



(58) and (59) follow from (20) and (^19) respectively. Since £2 ^5-2 < 1 for 13 > 5/6 and 
6 > e'^/\ is 0(1) in C{n) as well as in C{D{H^)) by (g. So (g follows from the 
fact that iff is uniformly bounded from D{Hq) to "H. Since 

and 6 > e^^^, (57) also follows from (65). Now (72) directly follows from Lemma [s] 

For (70) we apply exactly the same reasoning as for (71), however, with (21) instead 
of (20), which worsens the bound. □ 



Now Lemma [sj applied to P^^ with n = 1, 5 ~ 5i ~ £ + £2/3^2 and /?! ~ £2/^5; 

yields a projector P^c with all the properties claimed in Propositions |3] and |4j 
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4 Proofs of Lemmas 



Proof of Lemma [T] 

The statement about the potentials is standard. Using for examples the estimates 
contained in Proposition 1 and in the proof of Proposition 2 of jTenj we can show easily 
that 

■ Pj,y,ii)\\c{Do,H) = W'^iv^al) ■ Pj,yXi)\\c{Do,H) < C||^^SL, 

where 



and, given a function / G L^(]R^ C^), 

ll/IU := (11/1^1 ^^^||l2(M3^c2) + II/I|Z,2(K3^C2)) ^ • 

Using these explicit expressions we get then 

\\^^'\(^yj)-Pj,y,(4c{Do,n) < CAfi-' = C^^^^ = CAo. 

In the same way we have 

Ih A{ay,f :||^(^.^,^) < C\\vj^l\\l < CiA^^-'f = CAl 

Identical results hold for the coefficients of the Hamiltonian containing the nuclear co- 
ordinates, so all the coefficients in can be bounded with an e-independent bound in 
terms of Hq or H^^^. 

Proof of Lemma [2] 

Since Hf is nonnegative and since -ff^oi ® ^ 1 ® iff commute, we have that 

(ii^„i ® 1)" < ® 1 + 1 ® iff)" = {H^,r (73) 

for any n G N. To estimate tensor product operators in C{Dq), the following character- 
ization of this operator norm will be useful. 

Lemma 11. Let {H,D{H)) be self-adjoint and 

:= {ijj G D{H) I H''^ G D{H) for k = 1, . . .n - 1} 



be equipped with the graph norm 



j=0 
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Then {D"', \\ ■ \\d^) is a Banach space, 

:= {{H + I V' e ^} = ^" 

and for A e C{D'^,D"^) the operator norm ||^||£(d",_d'") is equivalent to the norm 

m 

i=o 

More precisely, there are constants Cn depending only on n ( not on H or A ), such that 



1 



l+m 



R{n,m) ^ ^ C'n||^||i?(n,m) 



Proof. Since H is self-adjoint, it is closed and therefore (-D", || ■ \\d") is a Banach space. 
Let ^ e D", then (i/+i)'^'0 e and thus ^ G L*^. Let conversely ^ = (i^+i)-"0 G D^, 
then if'^'^ eH ior k < n since if (if + i)~^ e C{T-i). For the norms observe that 

m m 
\\A^\\j,r. = J^llii^A^II = 5^||ii^'A(ii + i)-"(ii + i)>|| 

j=0 j=0 

m 

< J2\\H^A{H + ir-\\cin)\\{H + iri^\\ < CJA\\n^^,^)\mnn 

j=0 

and thus 

Conversely for j <m 

||ii^A(ii + i)->|| < ||A(ii + i)->||^^ < ||A|U(0„,0^)||(ii + i)->||0n 

< ||^||£(D",D-)||V'||, 

where we use ||(if + i)~^||£(DJ,DJ+i) = 1- Thus 

||^||ii(n,m) < ("^+ 1)||A||£(D",D-) ■ 

So the C{Dq, D^)-norm of an operator S (8) 1 for m < n is estimated by 

m 

\\B ® 1|U(D„",D,"') <CnJ2 IK^oT {B ® 1) + i)""ll ■ 

j=0 

This will turn out useful, since for bounded operators on Hilbert spaces ||A <S) B\\ — 



□ 
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• Let's look at a single term in the sum more closely, 

(H^.y (B ® 1) (H^, + i)-" = ® 1 + 1 ® H,y {B ® 1) {H^, + i)"" 

= E ((^moi)' ® 1) (l ® H]-') {B 1) (ifo^ + i)-" 



i=0 

j 



Since 



£=0 

X + iY^-^ ® (i/f + i)^-^) {H^, + i)-« . 

((i/^„i + i)^+"-^®(//f + irO(/7o' + i)-" and i/r'(^f + i)"^'"'^ 



are bounded uniformly in e due to (73), it suffices to control terms of the form 

By Lemma these are controlled again in terms of ||i?||^(.j-,<!+n-j y 
Proof of Lemma [3] 

The assumption ||A||£(Do,W) < 5 < 1 implies that foripED 

\m\<6\mno = smoi'\\ + \m) 

and thus A is i^o-bounded with relative bound smaller than 1. The equivalence of the 
norms follows from 

ii^iidh = ii^^^ii + ii^ii < iii^ov^ii + iiv^ii + \m\ < ii^iido(i + s) 

and 



The last claim follows from the Helffer-Sjostrand formula 

xiH) = - [ d,x{z){H~z)-Uz, 
where x is an appropriate almost-analytic extension of x, and the resolvent formula 

< 5{l + 5) (l 



\z\ + V ' 
|Imz| 
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Proof of Lemma 2] 

Due to the smearine; of the nuclear charge it holds Kn,14n e C{f (M^^, ^^(MSr^^^ ^^^^ 
that 

Thus the mapping x i— )■ {Hci{x) — z)~^ is in Cl^ (M^', C{l-Le\)) ■ Since Ej is separated by a 
gap, the projection Pj{x) associated to Ej{x) is given via the Riesz formula: 

Pj{x) = ^ (p (Helix) - z)~^ dz, 

where 7(x) is positively oriented closed curve encircling Ej{x) once. It can be chosen 
independent of x locally because the gap condition is uniform. Therefore {Hei{-) — z)~^ G 
Cl{R^\C{Hei)) entails that Pj G Cl{R^\ CiHei)) ■ By 

Ej{x)Pj{x) = Hei{x)Pj{x) = ^ (f z{H,y{x) - z)~^ dz 

J-y{x) 

we obtain EjPj G ^^(M^^ ^CHci))- Then Ej = ti L2(RSs){EjPj) G ^^^(M^O- For it holds 

V.tr(E,P,) = V.tr((E,P,)P,) = tr((V.E,P,)P, + {E,Pj)V.P,) 
= tT{{V,EjPj)Pj) + tr((E,P,)V,P,) < oo 

because Pj and EjPj are trace-class operators and the product of a trace-class operator 
and a bounded operator is again a trace-class operator (see e.g. |ReSii| , Theorem VI. 19). 
The argument for higher derivatives goes along the same lines. 

For the last claim we observe that if^(9"Pj- is bounded for any a G Ng' and any 
n G Nq. Since H^i, Ej and Pj have bounded and smooth derivatives, this can be easily 
seen inductively by differentiating the identity 



= iH,i~E,rP^. 



Proof of Lemma [5] 



We proceed by induction. For m = 1, i.e. \a\ < 2, we have by standard elliptic estimates 
that e'^'c?^ is relatively bounded by — e^A^,. Now H^^^ has the form 

^moi = -^^^x ® 1 + H,i{x) = -e^A.^ ® 1 + 1 ® (-Aj, + \/ec(y)) + Ven{x, y) , 

' V ' 

=://c,o>0 

where Vg^ is bounded with bounded derivatives. Hence 

||-e2A.^||<||i/^>|| + ||Kn||oo||^|| 

and thus 

\\A^e-d:^\\ < \\AJ \\e"d:^\\ < C{\\H^^,M\ + M) = 
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with a constant C independent of e. 

Now assume that we proved the assertion for operators of order n — 1 and let |a| = n 
and m = \n/2]. Then A{x)e^°'^d" is relatively bounded by (— e^Aa,)™ again by standard 
elliptic estimates. Using the induction hypothesis we find that 

||(5^A.r^|| < C{\\H^Ue'A^r-'^p\\ + \\{e'A^r-'^\\) 

since (e^A^)"*"-*^ and [H^^^, {e"^ Ax)"^~^] are both differential operators of order at most 
2m - 2 < n - 1. 

For the second claim note that 

since [{H^^^)'' , Aa]e"d" and [(H^^^)'', e"d"] are admissible of order 2/c — l + n < 2{k + m). 



Proof of Lemma [6] 

All the operators appearing are differential operators with coefficients Aa that are com- 
posed of derivatives of Pq and R, i.e. of O^Pq and d^R. So they are all admissible in 
the sense of Lemma HI It remains to show that also commutators of the coefficients 
Aa with (if^Qj)'^ are admissible, which in turn follows if H^d^Aa is bounded for any 
(3 e Nf. Now according to Lemma|4|9fPo(a;) e £{71^, D{H^^)) for any n and clearly also 
d^R{x) E C[D{H2)) for any n. Since every coefficient Aa appearing in the construction 
contains at least one factor of the type d^Po, the claim follows. 



Proof of Lemma [7] 



First take any G C^(]R) and an appropriate almost analytic extension. Then the 
Helffer-Sjostrand formula implies 

- [ d,^{z){H-z)-'[A,H]lj{H){H-z)-Uz 



< 



< 



C{'H,D(H")) 

\d,^iz)\\\iH - z)-'\\ciDiw.)) WiH - z)-'\\ciH) dz 

1 



|Im(2;) 



dz < CJ. 



Taking the adjoint shows that also \\li{H)[(f){H), A]\\c(n) < Ccf,^- With the bound 
\\'i-i{H)\\c(n,D{H")) < Cn we get also 

\\li{HmH),A]\\ccH,DiH-)) < CnC^S. 
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Now choose x ^ C'^f (M) with suppx C / and x\i = 1- This imphes x = X^i{H), 
X = x^i{.H) and xx = X- Using the above estimate for = x and (p = x, we get 

\\X{H)A-Axm\ = \\x{H)A-Ax{H)x{H)\\<\\x{H)A-x{H)Ax{H)\\+CJ 
Proof of Lemma [8] 

We first state another lemma that will be used in the proof. 

Lemma 12. Let {H,D{H)) be a self-adjoint operator and equip the domains := 
D{H'') with the graph norms \\^\\d^^ ■= ELo II^VII and let N e N. If A e C{H) 
satisfies 

\\[AH]\\c{D",D"--^) < Sn 

for all 1 < n < N , then 

n 

\\A\\c^D^^<n\\A\\ + J2^^ (74) 

i=l 

for all 1 < n < N . Moreover, if in addition 

\\{A-z)-^\\cin) < « 

and Sn < 1 22^' ^^s'^ 

\\iA-z)-'\\ciD-)<2'-a (75) 

for all n < N . 

Proof. We proceed by induction. 

71 n— 1 n—1 

^\\AHij\\Dr.-i + \\[H,A]i;\\nr.-i 



< Pll ll^ll + + \\[H,A]\\c{D"^D"-^)\\V\\D'^ 

and thus 

||^||/:(Z)") < + ll^ll + ■ 



Since for n = 1 the computation yields ||v4||£(£i) < ||y4|| + 6i, this implies (74). For (75) 
we proceed analogously and abbreviate R := {A — z)^^ . 

n n—1 n—1 

i=0 1=0 i=0 

= IWW + WRH^Wnr.-. + \\R[A, H]R^P\\nn-^ 

< \\R\\ WipW + ||i?||£(Dn-l)||^||£,n + 1 1 /2 1 1 1 ) 1 1 ^ 1 1 £(D" ) 1 1 1 1 



45 



and thus 



\R\\ciD") < 



l^ll + \\R\\C(D^ 



1 — 5n||-R||£(D"-l) 



For n = 1 the this yields ||-R||/:{i:)) < 4q; if l — 6ia > ^ and by induction one obtains (75). 

□ 



Sine Q is self-adjoint in C{'H), (50) implies that the spectrum of Q is located in balls 
of radius 6i around and 1. Thus, for 5 < ^ the curve 7 : [0, 2tt) — )■ C, 7(6') = 1 + |e^^, 
is contained in the resolvent set of Q G j^i'H) and we can define 



Q:- 



2tt 



as a bounded operator in C{l-L). Note that Q is just the spectral projection of Q related 



to the spectrum near 1. For simplicity write R{z\ 
for z G 7 we have < 4 and by Lemma 12 



{Q — z) ^ and assume 5 < \. Then 



for 5r,.< also 



2 22 



mm 



< 4 ■ 2 



2n 



(76) 



is uniformly bounded on 7. Hence < 4"+^ for all n < N. The fact that Q — Q 

has spectrum only in a ball around of size 6 implies 



\\Q-Q\ 



cm 



<6. 



To estimate the difference also in C{D^) we use Nenciu's formula |Nenj 



Q-Q 



i / R{z) - R{1 - z) 



2tt 



1-z 



dz (QQ-Q) . 



Now (50) together with (76) implies the second part of (52) with C„ = 8 ■ 4""*"^. From 



[i/,g] = — (1^ R{z)[H,Q]R{z)dz, 



and (51) it follows that 

II [if, Q] ||£(D",D"-l) < CnSn 

with Cn = 4^'^"'"^. From now on we will not keep track of the exact value of C„ and 
increase it as necessary in the following steps. But it should always be noted that it 
depends only on n. 

Now pick X G C^(]R) with support in {e — 1,E + 1) and with xl^+i = '^e+^- Then 
(53) together with Lemma [TI implies that there is a constant C depending only on x, 



which in turn can be fixed given E, such that 

\mH),Q]\\c(n,D-)<C(3^. 

Thus for 2; G 7 

\\[Riz),xm\\cin,D^) = II^W [xiH),Q] R{z) 
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<CnC(3^ 



which shows that 



[iJ,Q]V|(if) = ^j^R{z)[H,Q]R{z)dzx{H)l^^r{H) 
= ^ f^Riz)[H,Q]xiH)Riz)dzl^_,^iH) 
+ R{z) [H, Q] [Riz),xiH)] dz 1 (H) 



imphes 



[H,Q] 1^+1 



< CuCe Pi 

/;(w,D") 



Finally we get (55) using again Nenciu's formula and (54). 
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